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ABSTRACT 

Here we extend the exploration of significantly super-Chandrasekhar magnetised white 
dwarfs by numerically computing axisymmetric stationary equilibria of differentially 
rotating magnetised poly tropic compact stars in general relativity (GR), within the 
ideal magnetohydro dynamic regime. We use a general relativistic magnetohydro dy¬ 
namic (GRMHD) framework that describes rotating and magnetised axisymmetric 
white dwarfs, choosing appropriate rotation laws and magnetic field profiles (toroidal 
and poloidal). The numerical procedure for finding solutions in this framework uses 
the 3 + 1 formalism of numerical relativity, implemented in the open source XNS 
code. We construct equilibrium sequences by varying different physical quantities in 
turn, and highlight the plausible existence of super-Ghandrasekhar white dwarfs, with 
masses in the range of 2 — 3 solar mass, with central (deep interior) magnetic fields of 
the order of lO^^G and differential rotation with surface time periods of about 1 — 10 
seconds. We note that such white dwarfs are candidates for the progenitors of pecu¬ 
liar, over luminous type la super novae, to which observational evidence ascribes mass 
in the range 2.1 — 2.8 solar mass. We also present some interesting results related to 
the structure of such white dwarfs, especially the existence of polar hollows in special 
cases. 

Key words: stars: magnetic fields - white dwarfs - gravitation - MHD - stars: massive 
- super novae: general 


1 INTRODUCTION 

Recently, there has been a lot of interest in exploring 
the possible existence of massive compact objects, in par¬ 
ticular neutron stars and white dwarfs, based on the 
theory of stellar structure. This is primarily because of 
the observational evidences fo r massive neutron sta r bi¬ 
nary pulsars PSR J 1614-2230 (IPemorest et al. ]H) and 
PSR J0348+0432 (jAntoniadis et al. I l2013l j with masses 
1.97M0 and 2.QIM0 respect i vely, w here M© is the mass of 
the Sun. IWeissenborn et al~l (|2012l j examined the possibil¬ 
ity of very massive neutron stars in the presence of hyperons 
and th e conditions to obtain the same. IWhittenburv et ahl 
(I2OI2I I. based on quark-meson coupling model, showed that 
the maximum mass of neutron stars could be ~ 2Mq , when 
nuclear matter is in /^-equilibrium and hyperons must ap¬ 
pear. Apart from the exploration based on the equation 
of state (EoS), as mentioned above, neutron stars with 
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mass > 2Mo were shown to be possible by exploring ef- 
fects of magneti c fields, with centr al field Be ~ 10^® G 
dpiir et al. ] |2014 ISinha et al. 1l20l3h. and modification to 
Einstein’s gravity (lArapoglu et ahl 20141 : lAstashenok et al~l 
I2OI2I: ICheoim et, al. 112013) ! 


On the other front, the discovery of over-luminous 
type la supernovae (e.g. SN 2003fg, SN 2006gz, SN 2007if, 
SN 2009dc) argues for significantly super-Chandrasekhar 
progenitors, of mass 2.1 — 2.8M0, in such super- 
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the super-Chandrasekhar limiting mass of white dwarfs as 
a potential possibility. Mukhopadhyay and his collaborators 
initiated the exploration of the possible existence of such 
a super-Chandrasekhar limiting mass, first by exploiting 
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exploring the effects of mo difying Einstein’s gravity 
(|Das Mukhopadhyay I [ 201 5bl l. 

In fact, the existence of super-Chandrasekhar white 
dwarfs was shown to be possible many decades back. 
Ostriker and his collaborators attempted to model mag¬ 
netised and/or rotating white d warfs in the Newtonian 
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magnetised Newtonian white dwarf was reported then 
to be 2.26Mc:), with , how ever, equatorial radius 2510km 
(lOstriker Tassoul I Il969l b A couple of decades later, 

I Adam I (|l986l l obtained significantly super-Chandrasekhar 
Newtonian magnetised white dwarfs, but did not impose 
realistic constraints on the magnetic fields. Nevertheless, 
none of them could pinpoint whether there is a limit¬ 
ing mass or not, in the same spirit as the Chandrasekhar 
mass-limit for non-mag netised whi te dwarfs was proposed 
( Chandrasekhar I Il935l b However, IPas Mukhopadhyay I 
( 2013al i showed the existence of a new mass-limit for white 
dwarfs in the presence of magnetic field, but in the absence 
of rotation, with a simpler (constant) field configuration. In 
fact, after the initiation by the Mukhopadhyay-group, re¬ 
cently the topic of super-Chandrasekhar white dwarfs has 
been given renewed attention by the community. 

Nevertheless, apart from magnetic fields, rotational ef¬ 
fects also could increase the mass of white dwarfs above the 
Chandrasekhar-limit, as was sho wn bv lOstriker Hartwick I 

iBoshkavev et al.1 (|2013l b In fact, in general, a white 


dwarf could be rotating as well as magnetised. Its angular 
frequency, however, conventionally, is expected to be much 
smaller than that of a typical neutron star. Nevertheless, in 
the presence of high rotation, white dwarfs are expected to 
be oblate spheroids. On the other hand, in the presence of 
poloidally and toroidally dominated magnetic fields, the re¬ 
spective white dwarfs are expected to be oblate and prolate 
spheroids. Therefore, the overall shape of rotating magne¬ 
tised white dwarfs depends on the relative effects of rotation 
and the compo nents of magnetic field. _ 

Recently, IPas Mukhopadhyay I (|2015al i have ex¬ 
plored, using the (modified) XNS code, general relativis¬ 
tic magnetohydrodynamic (GRMHP) analyses of magne¬ 
tised but nonrotating white dwarfs and showed that their 
mass could be as high as 3.4M0. In this paper, we ex¬ 
plore GRMHP numerical modeling of white dwarfs, includ¬ 
ing both magnetic and rotational effects, by using the XNS 
code. Note that the XNS code is particularly suited to solv¬ 
ing for highly deformed st ars due to a strong magnetic held 
and/or rotation al effects (|Bucciantini X Pel Zann^ 1 20 111 : 
IPili et al~ll2014l i . Although originally the XNS code was de- 
veloped to investigate neutron s tars, recently it was modified 
(|Pas X Mukhopadhyay 112015al l in order to obtain equilib¬ 
rium configurations of deformed, magnetised white dwarfs. 
Here, we further modify it to make it suitable for deformed, 
magnetised, differentially rotating white dwarfs. 

This paper is organized as follows. In Section 2, we 
briefly mention the numerical approach as described by 
|Pas X Mukhopadhyay I (|2015al b In section 3, we describe 
the organisation of the results to be presented in the sub¬ 
sequent sections and chosen basic features/assumptions and 
units for the work. In section 4, we discuss the results of 


non-magnetised, differentially rotating white dwarfs. In sec¬ 
tions 5 and 6, we present the results for differentially ro¬ 
tating white dwarfs having toroidal and poloidal magnetic 
field configurations respectively. We also explore the possi¬ 
bility of magnetised, rotating white dwarfs having very small 
radii, in section 7. We finally conclude in section 8 with a 
summary. 


2 NUMERICAL APPROACH AND SET-UP 


We refer the rea ders to Bucciantini X Pel Zanna ( 201lh . 
IPili et af~l (|2014l l and |Pas X Mukhopadhyay I (|2015a ) for a 
complete description of the GRMHP equations and numer¬ 
ical approach. Our interest is obtaining equilibrium solu¬ 
tions of high density, rotating magnetised, relativistic white 
dwarfs. 

For the convenience, we recall below the equation per¬ 
taining to the differential rotation profile chosen in this 
paper. Stationary configurations can have many rotation 
laws consistent with the hydrostatic balance condition in 
GRM HP. One such form , deduced based on the Newtonian 
limit (|Stergioulas Il2003l b is used in XNS, given by 


j{Q) = (Qc - Q), (1) 

where j is related to the specific angular momentum, Q, is the 
angular velocity of a local zero angular momentum observer 
(ZAMO) as measured from infinity, Qc is the angular veloc¬ 
ity at the centre of the coordinate system, and A is a positive 
constant that is implicitly related to Qc- We will refer to A as 
the differential rotation parameter. When Q ^ Qc, A ^ 00 
and the system is of uniform rotation. In the opposite limit, 
when A ^ 0, j becomes constant in space. Finite values of 
A indicate a spatial variation of angular velocity, and hence 
differential rotation. 


3 ORGANIZATION OF RESULTS 

We have grouped the various sequences broadly under three 
sections — non-magnetised rotating sequences, sequences 
with toroidal magnetic field, and sequences with poloidal 
magnetic field. 

At the outset, we illustrate the two kinds of field geome¬ 
tries in Fig. [1] which will be considered throughout. It may 
be noticed that the toroidal held is zero at the centre and 
peaks at a small radius, distance r from the centre, within 
the star. All the chosen toroidal held prohles, throughout 
the paper, correspond to m = 1.4, unless stated otherwise, 
where the number m is the toroidal magnetisation index that 
determines t he form of the barotropic magn etic held pro- 
hle (see, e.g. jBucciantini X Pel Zanna Il201ll b The poloidal 
held we choose is purely dipolar, having no non-linear cur¬ 
rent sources. Various held prohles, chosen in the paper, are 
determined by our respective choice of the maximum value 
of magnetic held and the corresponding magnetic dipole mo- 
ment obtained for the self-consistency of the s olutions (see 

et al. I I 2 OI 4 [Pas X Mukhopadhyay Il2015al . for details). 

We compute all our sequences at the central density 
Pc = 1.9902 X 10^°g/cc, which is high enough to ensure al¬ 
most complete relativistic electron degeneracy, so that we 
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(a) (a) (b) (b) 


Figure 1. Illustration of magnetic field (B) geometries: The pan¬ 
els are contour plots of the magnetic surfaces, and B in units of 
for (a) purely toroidal field, and (b) purely poloidal field, 
when r is in the units of 1.48km. 



log(p) 

Figure 2. Polytropic P, that locally fits the Chandrasekhar EoS 
for a zero temperature electron gas, as a function of density in 
g/cc. 

may use the polytropic EoS with n = 3 (i.e. P = |) consis¬ 
tently throughout the star. This density is still lower than 
the limit at which gravitational instabilities in GR set in, 
which is about 3 x 10^°g/cc. 

The choice of a polytropic EoS with n = 3 throughout is 
almost appropriate for white dwarfs, as can be seen from Eig. 
[21 which shows the best poly tropic T, that fits approximately 
the full zero temperature EoS locally, at various densities. It 
is clear from the figure that T = 4/3 gives a good fit in a wide 
density range, from about 10^ to 10^^g/cc. Since the outer, 
low density, layers of a white dwarf do not contribute much 
to the global properties such as total mass, we can safely 
use this model to make conclusions about such properties of 
white dwarfs. 

Hence, all the sequences described in this paper have 
been computed with the fixed EoS 

P^kspK ( 2 ) 

where ks = 4.9018 x lO^^cc s~^ ■ 

In the discussion and all the tables and figures that 
follow, the units used are the following, unless mentioned 
explicitly. 

(i) Mass - solar mass (M©); M© = 1.987 x 10^^g. 

(ii) Radius - km. 

(iii) The radial coordinate in all the contour plots is in 
units of ro = 1.48km (chosen because the normalisation G = 
c = M© = 1 has been followed). 

(iv) Magnetic fields {B) - lO^^G. 

(v) Angular velocities (Q) - rad/s. 

(vi) The normalisation factor of density for the contour 
plots of log(p/po) is po = 10^°g/cc. 


All the chosen differential rotati ons (jKomatsu et al. |[l989al : 
iBucciantini Sz Del Zanna1l201lh correspond to — 10^, 
unless stated otherwise. 

Numerical solution of equations involves discretisation 
of the domain. We have used a resolution of Nr = 250 for 
the radial coordinate, and Ne = 100 for the angular coordi¬ 
nate. The XNS solver uses a compactified domain, so that 
the range for the radial coordinate always starts at 0 and 
truncates at some specified Rmax- Eor the present sequences, 
Rmax varies from 1500km to 4000km. It may be noted that 
for numerical purposes, the pressure and density at the sur¬ 
face of the star are not zero, but a small, finite value at 
which the integration is truncated. 

The choice of gridjWe have followed, is drawn from the 
guide to the XNS codelj The angular coordinate varies from 
^ = 0 to TT. The choice of 100 angular grid points appears 
to be sufficient to capture the variation of quantities with 
0 for the models presented here, which have no cusps or 
other sharp features. The value of Rmax can, in principle, 
be arbitrarily chosen. However, it is necessary to guarantee 
that the white dwarf is well resolved over a sufficient number 
of grid points (50-100), so this parameter and Nr need be 
chosen consistently. Our choice of Nr = 250 ensures that all 
the configurations we compute have well over 150 grid points 
within the stellar surface. The Tolman-Oppenheimer-Volkoff 
(TOV) equation solver routine of XNS, which computes the 
initial guess TOV solution for the metric solver, converges 
when the ADM masses measured at Rmax and Rmax I co¬ 
incide within a given tolerance. Eurthermore, the iterative 
metric solver of XNS, based on conservative variables, mea¬ 
sures convergence by the quantity D = which contains 

the unknown value of the conformal factor -0 at the centre. 
All the results presented in this paper correspond to conver¬ 
gent models and are stable at the selected resolution. 

We also recall at the outset that the non-magnetised, 
non-rotating spherical equilibrium configuration of a 
white dwarf with the central density pc = 1.9902 x 

10^°g/cc has mass M — 1.4 16 and radius Ve = 1215 
(|Das V Mukhopadhvav Il2015al i. 


4 NON-MAGNETISED ROTATING 

CONFIGURATIONS 

We start with the sequences of non-magnetised rotating 
stars, focussing on differentially rotating configurations. 
Uniform rotation has a marginal effect on the mass and ra¬ 
dius of the sta r, except close to the mas s-shedding limit. As 
pointed out bv lOstriker X Mark I (|l968h . uniformly rotating 
sequences are terminated at low values of the kinetic to grav¬ 
itational energy ratio (KE/GE=0.05-0.06), by the balancing 
of gravitational and centrifugal accelerations at the surface. 
However, uniform rotation does result in severe distortion of 
the outer laye rs. _ 

lOstriker! (|l970h arrived at a criteria for stability of ro¬ 
tating stars based on the value of KE/GE, by showing that 
configurations become unstabl e to non-axisymmetric p ertur- 
bations when KE/GE ^ 0.27. iKomatsu et al. I (|l989al i have 


^ http://www.arcetri.astro.it/science/ahead/XNS/files/Guide.pdf 
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adopted a similar criterion for stability of rotating configura¬ 
tions, and argued that when relativistic effects are stronger 
and ring-like or toroidal structures form, the KE/GE thresh¬ 
old reduces to 0.14 — 0.16. With this in mind, we have 
always terminated our sequences at KE/GE ^ 0.14, even 
when the code is capable of converging to an equilibrium 
solution for higher values of this ratio. 

We plan to explore the phenomenon of ‘polar hol¬ 
lows’. The term polar hollow refers to a change of the den¬ 
sity isocontours of the star from a convex to a concave 
nature, near its poles along the rotation axis. Polar hol- 
lows in polytropic stars have previously been reported by 
iKomatsu et al. [ (Il989blb and investi g ated in detail for the 
Newtonian case bv ISmith Gollinsi (|l992l b To the best of 
our knowledge, not much work has been done on polar hol¬ 
lows in rotating and magnetised stars in GR. In this work, 
we take a first step towards understanding this problem, 
by computing and demonstrating the plausible existence of 
equilibrium solutions with polar concavities, for the specific 
angular momentum distribution profile given by equation 
©• 


4.1 Differentially rotating sequence with varying 

Eixing the differential rotation parameter A (see §2.2 of 
iKomatsu et al. Nl989al . for details of the physical meaning 
of A), we vary the central angular velocity Qc- The stel¬ 
lar density profiles in Eig. [3] show the development of polar 
hollows with the increase in Qc- In fact, the polar radius 
(vp) decreases, whereas the equatorial radius (re) increases, 
with the increase of Qc, as might be expected due to in¬ 
creased centrifugal forces. This results in oblate structures. 
The appearance of polar hollow is due to the differential 
rotation of the star. A faster rotation in the inner region 
disperses the matter to the outer region, making the central 
region of the star very oblate. However, the white dwarfs 
become less oblate in the outer region, with the decrease of 
rotational frequency, hence develops a polar hollow. Eor a 
KE/GE=0.14, and equatorial angular velocity at the outer 
surface Qeg = 3.089 corresponding to a rotational time pe¬ 
riod ~ 2s, the mass is about 1.826 — an increase of nearly 
29% from the non-rotating case mass of 1.416. Table 1 shows 
the effects of increasing Qc and Qeq in the white dwarfs. 

In Eig. m we show the equilibrium sequences of vari¬ 
ous stellar quantities. There appears to be no natural limit, 
such as the Keplerian limit for uniform rotation, for differ¬ 
entially rotating objects with the rotation profile of equa¬ 
tion ©• The ratio KE/GE apparently increases without 
bound. However, instabilities are expected to terminate the 
sequence at KE/GE 0.27. 

Eor comparison, we note that the uniformly rotating 
configuration of highest Q, that we were able to compute, 
has the features: Qc = 5.778, M = 1.455, re = 1766, Tp/ve = 
0.672 and KE/GE=0.0176. This is actually quite close to the 
mass-shedding limit. In fact, naively computing 
where ge = GMjrl, we find it to be ~ 0.95. We were unable 
to bring XNS to converge to a solution at higher values of 
Qc at the central density chosen here. 



Figure 3. Sequence of non-magnetised differentially rotating con¬ 
figurations with changing Qc- The panels are contour plots of 
log^^^ corresponding to the Qc values (a) 2.028, (b) 12.168, (c) 
22.308, (d) 32.448, for which various physical quantities are listed 
in Table 1. 






(c) (d) 



(e) 


Figure 4. Non-magnetised differentially rotating sequence with 
changing He. The panels depict (a) mass, (b) equatorial radius, (c) 
KE/GE, (d) equatorial angular velocity, (e) r^/re, as functions 
of central rotation rate. 
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Qc 

M 

Te 

Qeg 

KE/GE 

rp/re 

2.028 

1.417 

1223 

0.257 

6 X 10-4 

1.000 

12.168 

1.465 

1258 

1.466 

0.022 

0.915 

22.308 

1.593 

1341 

2.401 

0.072 

0.753 

32.448 

1.826 

1435 

3.089 

0.140 

0.588 


Table 1. Non-magnetised differentially rotating configurations 
with changing Qc- 


A2 

M 

Te 

Qeg 

KE/GE 

TpIVe 

0 

1.414 

1217 

0 

0 

1.000 

5 X 10® 

1.438 

1288 

2.2838 

0.0115 

0.927 

2.5 X 10® 

1.448 

1406 

4.234 

0.0158 

0.832 

1- 

o 

1—1 

1.451 

1524 

5.219 

0.0171 

0.767 


Table 2. Non-magnetised differentially rotating configurations 
with changing r2c = 5.778 fixed. 



4.2 Differentially rotating sequence with varying 

We fix Qc = 5.778, which is close to the mass-shedding limit 
for the uniformly rotating configuration of the chosen central 
density. Figure [5] and Table 2 show that as ^ oo, we 
recover the uniformly rotating configuration (see ^4.ip . We 
note that the sequences of stars in Figs. [5] and [6] show no 
sign of polar hollows. Our studies indicate that polar hollows 
develop only when Qc is larger than the angular velocity of 
the corresponding uniformly rotating model at the mass- 
shedding limit. 


5 MAGNETISED CONFIGURATIONS WITH 
TOROIDAL FIELD 

We choose a base configuration, by looking for a case where 
the ratio of magnetic to gravitational energies (ME/GE) 
and KE/GE are nearly equal, around which to construct 
the various sequences. The motivation behind this is that 
such a situation is roughly indicative of an equipartition of 
the effects on the white dwarf structure between the rota¬ 
tion and the magnetic field. Our base configuration has the 
features: maximum magnetic field Bmax = 3.227, m = 1.4, 
Qc = 30.42, Qeg = 0.988, ME/GE=0.144, KE/GE=0.136, 
M = 2.586, and Ve = 2543. 

Stability criteria for magnetised and rotating stars in 
GR are not well known. This is an important question that 
we plan to take up in a later work. How ever, based on the 

S iption for twisted-torus geometries, iGiolfi &: Rezzollal 
showed that compact stars with highly toroidally 
dominated magnetic field could be stable and realistic. Eor 
the sequences that follow, we choose to keep ME/GE below 
0.2, purely as a tentative criterion. It is, however, known 
that when the magnetic field exceeds 5, the EoS is 
strictly modified due to the eff ects of Landau quantisation 
(jPas Mukhopadhyav ir2012al L and hence we do not allow 
the field to go beyond 5 in magnitude, as we plan to stick 
with F = 4/3 throughout. 


Figure 5. Sequence of non-magnetised differentially rotating con¬ 
figurations with changing A^ and Qc = 5.778 fixed. The panels 
are contour plots of log^^^ corresponding to the A^ values (a) 

0, (b) 5 X 10^, (c) 2.5 X 10®, (d) 10^, for which various physical 
quantities are listed in the Table 2. 



log(A5 

(a) 



(b) 



log(A5 



0 1 2 3 4 5 6 7 

log(A5 


(c) (d) 



log(A5 


(e) 


Figure 6. Non-magnetised differentially rotating sequence with 
changing A^, at fixed Qc = 5.778. The panels depict (a) mass, (b) 
equatorial radius, (c) KE/GE, (d) equatorial angular velocity, (e) 
rp/re, as functions of the rotation parameter A^. 
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Bmax 

M 

Be 

ME/GE 

Bp! Be 

0 

1.416 

1215 

0 

1 

2.028 

1.510 

1358 

0.038 

1.010 

2.822 

1.715 

1639 

0.098 

1.032 

3.405 

2.238 

2330 

0.194 

1.068 


Table 3. Non-rotating configurations with purely toroidal mag¬ 
netic field, with changing Bmax- 



-1000 0 1000 



0 



-1000 0 1000 



0 


5.1 No rotation 

The physics of magnetised white dwarfs is complicated by 
the fact that the pressure becomes anisotropic due to the 
inclu sion of a magnetic field pressu re and a magnetic ten¬ 
sion (|Das Mukhopadhyay I [ 20 14bl b The magnetic tension 
for a toroidal field tends to cause the star to become prolate, 
stretching it along the polar direction. This is indeed seen 
from the density isocontours in Fig. [3 and Table 3. We also 
observe that the magnetic field does not influence the struc¬ 
ture of the star significantly until it attains a magnitude of 
about 10^^ G. This can be understood physically from the 
fact that the magnetic pressure term, becomes com¬ 

parable to the matter pressure in a white dwarf of central 
density ~ 10^°g/cc, with the EoS given by equation at 
this magnitude of field only. 

Figure [8] shows that with increasing field, both Ve and 
Vp increase, but Vp increases more than Ve leading to prolate¬ 
ness. However, the apparent structure is still nearly spher¬ 
ical, as Bp/re is atmost 1.07. We see that for a Bmax of 
about 3.4, the mass increases to 2.238, by about 58% from 
1.416 in the non-magnetised case. However, the value of re, 
2330, is nearly double that in the non-magnetised case, so 
that the mean density in comparison goes down by a fac¬ 
tor of nearly 10. This also aligns well with the fact that the 
density isocontours show a decrease in the degree of central 
condensation. 

Our computations reveal that ME/GE seems to increase 
without bound, with the code still capable of converging to 
an equilibrium solution for values of ME/GE even higher 
than 0.3, for example. This makes the need for a stability 
analysis all the more significant, and we also plan to tackle 
this problem in a future work. 


5.2 Uniform rotation 

We seek to construct a sequence of varying Qc, keeping the 
magnetic field fixed. In practice, although we fix the two 
constants that determine the magnetic field profile, there is 
a small change in the field with the change of Qc- This is due 
to the computation of a self-consistent solution. The effect is 
very small, the fractional change being ~ 10“^. The initial 
fixed Bmax is 3.092, and becomes 3.094 when Qc is increased 
from 0 to 3.4. 

Only in the far outer layers, where the centrifugal term 
(proportional to r) becomes large, the density isocontours 
become oblate — the deformation typical of axisymmetric 
rotation. The interior layers show mainly a magnetic field 
induced deformation, which is prolate for toroidal fields, as 
shown in Fig.[9l ME/GE remains nearly constant at 0.133 as 



Figure 7. Sequence of non-rotating configurations with a purely 
toroidal magnetic field, with changing Bmax - The panels are con¬ 
tour plots of logcorresponding to the Bmax values (a) 0, (b) 
2.028, (c) 2.822, (d) 3.405. The corresponding physical quantities 
are listed in Table 3. 




(a) 


(b) 



1.07 

1.06 

1.05 

1.04 

1.03 

1.02 

1.01 

1 


iog(B^ 


(c) 


(d) 


Figure 8. Non-rotating sequence with changing Bmax of toroidal 
magnetic field. The panels depict (a) mass, (b) equatorial radius, 
(c) ME/GE, (d) rpjre, as functions of the maximum magnetic 
field. 


Qc is changed, and KE/GE increases from 0 to about 0.015, 
as shown in Fig. [TOl 

The ratio QlqVelge comes out to be ~ 0.93 for Qc = 
3.407. The saturation of the effect of uniform rotation can 
be seen in the drastic decrease of r^/re, from 1.038 to 0.696, 
approaching the value at the mass-shedding limit for the 
non-magnetised case, as given in Table 4 and shown in Fig. 
M The mass increase due to the effect of rotation on the 
initial non-rotating magnetised case is about 0.05, which ap¬ 
pears to be the same as for the non-magnetised case, given 
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Qc 

M 

re 

KE/GE 

ME/GE 

rp/re 

0 

1.878 

1869 

0 

0.1342 

1.038 

1.358 

1.885 

1905 

0.002 

0.1339 

1.019 

2.620 

1.905 

2100 

0.008 

0.1333 

0.916 

3.407 

1.928 

2738 

0.015 

0.1326 

0.696 


Table 4. Uniformly rotating configurations with purely toroidal 
magnetic field, with changing Qc and Bmax = 3.092 fixed. 
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(b) 
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Figure 10. Uniformly rotating sequence with changing Qc at 
fixed toroidal Bmax = 3.092. The panels depict (a) mass, (b) 
equatorial radius, (c) ME/GE, (d) KE/GE, (e) r^jrc^ as functions 
of the central rotation rate. 


5.S.1 Sequence with fixed Bmax constructed by varying Qc 


Figure 9. Sequence of uniformly rotating configurations with a 
purely toroidal magnetic field with changing Qc at a fixed mag¬ 
netic field having Bmax = 3.092. The panels are contour plots 
of log^^^ corresponding to the values (a) 0, (b) 1.358, (c) 
2.620, (d) 3.407. The corresponding physical quantities are listed 
in Table 4. 

by Table 2. However, note that due to the problem of resolu¬ 
tion, Figs. fTOT b) and (e) do not appear to be smooth (show 
discrete behaviour, which is not physical). At higher resolu¬ 
tions, much larger runtimes are taken to compute a configu¬ 
ration, without revealing any meaningful physical change in 
the results. However for the primary understanding of the 
problem, our resolution suffices and the trends obtained are 
clear. 


5.3 Differential rotation 

Physically, differential rotation interacts with magnetic 
fields through rotationally induced currents. Hence, there 
is a rich range of phenomena associated with rotating mag¬ 
netised systems, such as magn etic braking and relati vistic 
jets from accreting systems etc. IMoss Smith ]^M) give 
a wide overview of the work that has been done on these 
effects. 


As before, the magnetic profile is fixed with m — 1.4, but 
the magnitude of Bmax slightly changes with Qc. As Ta¬ 
ble 5 and Fig. [11] show, for the range of results that we 
present based on our tentative stability criteria, Bmax does 
not change much. However, the Fig. [12] shows that Bmax, 

re, KE/GE, ME/GE, ah have quite a sharp increasing 
trend near the highest value of Qc = 50 in our sequence. 
This might be indicative of some limiting behaviour, and 
needs to be investigated analytically. 

We note from Eig.[TT]that the central regions of the stars 
are prolate, with strong polar hollows at higher Qc. Along 
the sequence, rp/re decreases, falling sharply near Qc = 50. 
The outer regions are relatively oblate with increasing Qc. 
The polar hollow appears to be sharper as compared to a 
case with no magnetic held but similar (but slightly smaller) 
rotation rate. The toroidal held peaks slightly away from 
the centre of the star, and the distortion of the outer den¬ 
sity isocontours to a prolate nature is pronounced at this 
radius. This could explain why the hollows, which are in¬ 
duced primarily by differential rotation, are accentuated by 
the toroidal held. 

Interestingly, Qeq increases at hrst, peaks and starts de¬ 
creasing. The mass increases by ~ 32%, from 1.878 for the 
non-rotating conhguration to 2.478 for a KE/GE value of 
0.121, and equatorial rotation period of about 6.5 seconds, 
as shown in Table 5. The equatorial radius increases to 2454, 
marginally more than the value for a held of similar magni¬ 
tude but without rotation. 
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M 

re 

^eq 

Bmax 

KE/GE 

ME/GE 

Vp/Ve 

0.003 

1.878 

1869 

0.0002 

3.0921 

1.5x10-9 

0.134 

1.038 

8.112 

1.918 

1922 

0.458 

3.1016 

0.011 

0.135 

0.982 

18.252 

2.097 

2118 

0.843 

3.1407 

0.054 

0.137 

0.816 

28.392 

2.478 

2454 

0.988 

3.2098 

0.121 

0.143 

0.617 


Table 5. Differentially rotating configurations with purely 
toroidal magnetic field with changing r2c and Bmax ~ 3.1 fixed. 



Figure 11. Sequence of differentially rotating configurations with 
a purely toroidal magnetic field with changing Qc and Bmax ~ 3.1 
fixed. The panels are contour plots of log^^^ corresponding to 
the Qc values (a) 0.003, (b) 8.112, (c) 18.252, (d) 28.392. The 
corresponding physical quantities are listed in Table 5. 

5.3.2 Sequence with fixed Qc constructed by varying Bmax 

For a Bmax of 3.584, the mass increases to about 3.159, by 
nearly 78% from the mass in the non-magnetised case. The 
radius behaves in nearly the same way as for the purely mag¬ 
netised case, and re becomes about 3322, more than double 
that of the non-magnet ised configuration. The surface Qeq 
reduces from 2.99 to 0.593, as Bmax increases. This could be 
understood as the increase of re increases centrifugal forces 
at the equator for a given Qeq, which enforces decreased 
Qeq in order to obtain equilibrium solutions. See Table 6 for 
details. 

We note again that the polar concavities, which develop 
primarily due to differential rotation, are accentuated by 
the toroidal field. Figure [13] clearly shows the severe polar 
denting with the increase in Bmax- 

An interesting feature seen in Fig. [THis that with the 
increase of Bmax^ KE/GE increases initially, attains a peak 
value of 0.136, and then starts decreasing. The value of r^/re 
too follows the corresponding pattern with a dip, when the 
largest KE/GE corresponds to the smallest r^/re- We spec¬ 
ulate that this may indicate the onset of some instability, 
and a possible upper limit to the mass in this sequence as 
well, and plan to follow it up in a future work. 



Figure 12. Differentially rotating sequence with changing Qc 
at fixed toroidal Bmax ~ 3.1. The panels depict (a) mass, (b) 
equatorial radius, (c) ME/GE, (d) KE/GE, (e) equatorial angular 
velocity, (f) r^/re, (g) Bmax, as functions of the central rotation 
rate. 


Bmax 

M 

re 

Qeq 

KE/GE 

ME/GE 

TpIVe 

0 

1.769 

1410 

2.990 

0.126 

0 

0.613 

2.299 

1.959 

1676 

2.180 

0.132 

0.046 

0.603 

2.996 

2.318 

2171 

1.339 

0.136 

0.108 

0.583 

3.584 

3.159 

3322 

0.593 

0.132 

0.203 

0.584 


Table 6. Differentially rotating configurations with purely 
toroidal magnetic field, with changing Bmax and Qc = 30.42 
fixed. 

6 MAGNETISED CONFIGURATIONS WITH 
POLOIDAL FIELD 


lEerraro I (|l937f ) showed that for a steady state, in a rotating 
star having a poloidal magnetic field, each field line must 
rotate uniformly. This result, called the ‘isorotation law’, 
dictates that the most general velocity field for a rotating 
star can only have uniform rotation of each poloidal flux 
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Figure 13. Sequence of differentially rotating configurations with 
a purely toroidal magnetic field, with changing Bmax and fixed 
r2c = 30.42. The panels are contour plots of log^^^ correspond¬ 
ing to the Bmax values (a) 0, (b) 2.299, (c) 2.996, (d) 3.584. The 
corresponding physical quantities are listed in Table 6. 
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Figure 14. Differentially rotating sequence with changing 
toroidal Bmax and fixed Clc = 30.42. The panels depict (a) mass, 
(b) equatorial radius, (c) ME/GE, (d) KE/GE, (e) equatorial an¬ 
gular velocity, (f) r^/re, as functions of the maximum magnetic 
field. 


Bmax 

M 

re 

ME/GE 

rpjre 

0 

1.416 

1215 

0 

1 

0.846 

1.424 

1210 

0.0046 

0.99 

2.260 

1.476 

1163 

0.0264 

0.91 

3.942 

1.610 

1045 

0.0721 

0.79 


Table 7. Non-rotating configurations with purely poloidal mag¬ 
netic field with changing Bmax- 

loop. This is a strong restriction, and makes the computation 
of models with poloidal fields more difficult. 

Poloidal fields are seen to result in highly centrally con¬ 
densed and compact configurations. Both re and rp tend to 
decrease with the increase of Bmax^ with the development 
of oblateness (equatorial bulge) due to a faster decrease of 
rp than re- In this respect, poloidal magnetic fields have an 
effect similar to rotation. 

Our sequences having a purely poloidal magnetic field 
have been constructed around the base configuration (rea¬ 
sons for this choice are discussed in the beginning of §5) 
with the features: Bmax — 4.575, Qc — 30.42, Vteq — 4.950, 
ME/GE=0.0773, KE/GE=0.106, M = 1.983, and re = 
1054. 


6.1 No rotation 

Eigures El and [16] show how the stars become progressively 
more oblate with increasing Bmax- ME/GE is lower than 
that for a toroidal field of similar magnitude. The mass in¬ 
creases from 1.416, in the non-magnetised case, to about 
1.610 with a Bmax of 3.942 — an increase of ~ 14%, as 
given in Table 7. However, like Eigs. llOT b) and (e), here also 
Eigs. Elb) and (d) exhibit the resolution problem (show¬ 
ing discrete behaviour, which is not physical), but clearly 
illustrate the trend of the quantity along the sequence. 

It is in general observed that the poloidal field con¬ 
tributes much less to mass increase than a toroidal field. 
However, it causes decrease in the size of the star - re de¬ 
creases from 1215, in the non-magnetised case, to about 1000 
for Bmax = 4. This results in more centrally condensed (hav¬ 
ing higher mean density) and compact configurations. The 
polar radius also decreases rapidly, causing oblate structures 
with rp/re as small as 0.79. 

The contrasting effects of poloidal and toroidal fields 
- poloidal causing greater central condensation and smaller 
radii, and toroidal causing reduced mean densities and larger 
radii - are entirely due to the difference in the field ge¬ 
ometries. The geometry decides the d irection and action o f 
Lorentz forces due the magnetic field (jGardall et al. llioOll l. 


6.2 Uniform rotation 

We notice in this case from Table 8 that including rotation 
in the model with a poloidal field seems to cause an abrupt 
decrease of mass initially, to 1.5, from 1.58. At higher rota¬ 
tion rates, it then rises marginally above the mass for the 
non-rotating star, to 1.62. However, the sequence numeri¬ 
cally terminates at similar ME/GE and mass as the mag- 
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M 

Te 

KE/GE 

ME/GE 

TpIVe 

0 

1.582 

1073 

0 

0.064 

0.803 

2.432 

1.509 

1085 

0.003 

0.057 

0.806 

5.159 

1.541 

1085 

0.012 

0.059 

0.749 

8.518 

1.625 

1420 

0.035 

0.064 

0.563 


Table 8. Uniformly rotating configurations with a purely poloidal 
magnetic field with changing and Bmax = 3.42 fixed. 



Figure 15. Sequence of non-rotating configurations with a purely 
poloidal magnetic field with changing Bmax- The panels are con¬ 
tour plots of log^^^ corresponding to the Bmax values (a) 0, (b) 
0.846, (c) 2.260, (d) 3.942. The corresponding physical quantities 
are listed in Table 7. 



Figure 16. Non-rotating sequence with changing Bmax of 
poloidal magnetic field. The panels depict (a) mass, (b) equato¬ 
rial radius, (c) ME/GE, (d) rpjre: as functions of the maximum 
magnetic field. 



Figure 17. Sequence of uniformly rotating configurations with 
a purely poloidal magnetic field, with changing Clc at a fixed 
magnetic field having Bmax = 3.42. The panels are contour plots 
of log^^^ corresponding to the r2c values (a) 0, (b) 2.432, (c) 
5.159, (d) 8.518. The corresponding physical quantities are listed 
in Table 8. 

onds, as shown in Figs. [17] and [181 This must mean that 
rotation has a stronger effect than the magnetic field, since 
the latter tends to cause a reduction of radius. The radius 
ratio decreases much more than its value of 0.67 (not shown 
in any sequence) and 0.79, in the case of pure rotation (not 
shown in any sequence) and pure magnetic field respectively, 
to 0.56. This is because the poloidal field and rotation both 
reinforce each other in making the star oblate. 


netised non-rotating configuration, so that we are unable to 
pursue the change of Qc furthermore. The magnitude of the 
field increases from 3.42 to 3.65, even though the param¬ 
eters determining the field are held fixed, indicating that 
the self-consistent solution requires non-negligible changes 
in the magnetic field for small changes (or no change) in the 
rotation profile. 

Moreover, the Ve increases, changing from 1073 in the 
non-rotating case to 1420, for a rotation period of 0.74 sec- 


6.3 Differential rotation 

In general, the allowed Qeg for rotating configurations with 
poloidal fields is much higher than in the similar cases with 
toroidal fields. This may be due to the decreased re in the 
former cases, which allows for faster rotation at the equator, 
because centrifugal effects decrease when radius decreases 
for a given angular velocity. 
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Figure 18. Uniformly rotating sequence with changing Qc at 
fixed poloidal Bmax = 3.42. The panels depict (a) mass, (b) equa¬ 
torial radius, (c) ME/GE, (d) KE/GE, (e) rpjre, as functions of 
the central rotation rate. 


6.3.1 Sequence with fixed poloidal Bmax constructed by 
varying Qc 


We note first that for self-consistency, Bmax, fixed initially, 
again increases from a value of 3.4, in the non-rotating con¬ 
figuration, to 4.8, for an equatorial rotation period of about 
1.15 seconds, as shown in Fig. [20] The radius changes little, 
due to a cancellation of the effects of rotation and poloidal 
magnetic field - the tendency of rotation to increase re, and 
that of the poloidal magnetic field to decrease the radius, al¬ 
most exactly balance each other in the range of parameters 
that we study. The mass increases by ~ 37%, from 1.502 to 
2.054, as shown by Table 9. The mean density increases by 
a factor of 2.3 from the non-rotating, non-magnetised case, 
reflecting the drastic increase in central condensation, and 
decrease in size. The formation of polar hollows is appar¬ 
ently hindered by the poloidal magnetic field, shown in Fig. 

M 

Figure [20] furthermore shows that the equatorial rota¬ 
tion rate increases as Qc increases — larger rates being al¬ 
lowed due to the smaller radius. ME/GE increases from 0.06 
to 0.08 along the sequence. However, at the higher values of 
Qc, KE/GE overtakes ME/GE, suggesting rotationally dom¬ 
inated effects. This is in contrast to the observed decreasing 
trend in radius, and the increasing central condensation. We 
speculate that this may be a purely GR effect, analogous to 
the case of a Kerr black hole, where increased spin results 
in more compact event horizons. Since KE/GE is greater 
than ME/GE at high masses, the mass increase must be 
more due to rotation than magnetic field effects. Neverthe- 



M 

re 


Bmax 

KE/GE 

ME/GE 

Tp/re 

2.028 

1.502 

1072 

0.322 

3.419 

6 X 10“^ 

0.057 

0.818 

12.168 

1.568 

1072 

1.929 

3.574 

0.020 

0.06 

0.769 

24.336 

1.798 

1072 

3.854 

4.119 

0.074 

0.071 

0.636 

32.448 

2.054 

1037 

5.429 

4.755 

0.117 

0.080 

0.556 


Table 9. Differentially rotating configurations with purely 
poloidal magnetic field with changing Qc and initial Bmax =3.1 
fixed. 



Figure 19. Sequence of differentially rotating configurations with 
a purely poloidal magnetic field with changing Qc and Bmax = 
3.1 fixed. The panels are contour plots of log^^^ corresponding 
to the Qc values (a) 2.028, (b) 12.168, (c) 24.336, (d) 32.448. The 
corresponding physical quantities are listed in Table 9. 

less. Fig. [20f F) reveals some unphysical resolution problem 
(with discrete behaviour). 


6.3.2 Sequence with fixed Qc constructed by varying 
poloidal Bmax 

Figure [21] and Table 10 show how the increase in poloidal 
Bmax results in a decrease of both re and r^, counteract¬ 
ing the effect of rotation. The polar hollows clearly become 
smaller, and we speculate that they may be entirely absent 
for high enough field strengths (which we were unable to 
compute due to convergence problems), if other instabili¬ 
ties do not disrupt the configuration. Figure [2^ (with some 
issues of resolution in Figs. [22] (b) and (e), as mentioned ear¬ 
lier) shows that KE/GE decreases from 0.125 to 0.103 with 
the increase of Bmax, indicating the restriction on rotation 
due to the virial theorem (energy conservation) and, hence, 
isorotation law. The mass increases by ~ 19.5%, from 1.695 
to 2.025, for a field of Bmax = 4.967, much less than, how¬ 
ever, the increase for a toroidal field of similar magnitude. 
The radius ratio decreases to 0.565, reflecting the extremely 
oblate configurations in Fig. 1211 The mean density rises by 
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Figure 20. Differentially rotating sequence with changing Clc and 
fixed initial poloidal Bmax =3.1. The panels depict (a) mass, (b) 
equatorial radius, (c) ME/GE, (d) KE/GE, (e) equatorial angular 
velocity, (f) r^/re, (g) Bmax, as functions of the central rotation 
rate. 


Bmax 

M 

re 

Heg 

KE/GE 

ME/GE 

TpIVe 

0 

1.695 

1409 

2.996 

0.125 

0 

0.610 

0.992 

1.710 

1373 

3.135 

0.124 

0.005 

0.613 

2.717 

1.805 

1231 

3.800 

0.118 

0.033 

0.597 

4.967 

2.025 

1019 

5.238 

0.103 

0.087 

0.565 


Table 10. Differentially rotating configurations with purely 
poloidal magnetic field, with changing Bmax and He = 30.42 
fixed. 

a factor of 2.4, showing the increased degree of central con¬ 
densation. 



Figure 21. Sequence of differentially rotating configurations 
with a purely poloidal magnetic field with changing Bmax and 
He = 30.42 fixed. The panels are contour plots of log^^^ corre¬ 
sponding to the Bmax values (a) 0, (b) 0.992, (c) 2.717, (d) 4.967. 
The corresponding physical quantities are listed in Table 10. 
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Figure 22. Differentially rotating sequence with changing 
poloidal Bmax and He = 30.42 fixed. The panels depict (a) mass, 
(b) equatorial radius, (c) ME/GE, (d) KE/GE, (e) equatorial an¬ 
gular velocity, (f) r^/re, as functions of the maximum magnetic 
field. 
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Figure 23. Illustration of smaller white dwarfs: The panels are 
contour plots of log for differentially rotating white dwarfs 

with poloidal magnetic fields having (a) M = 1.95, Ve = 746, 
Tp/re = 0.564, and (b) M = 1.72, re = 428, rp/re = 0.689. For 
other physical quantities, see the description in the text. 

7 EXPLORATION OF WHITE DWARFS WITH 

SMALL RADII 

As discussed above, white dwarfs having poloidal magnetic 
fields can have re and rp much smaller than those having 
toroidal magnetic fields. The question is, how small could a 
white dwarf be in size? In order to explore the answer to 
such a question, we consider the following two cases. 

(a) A white dwarf with pc = 3 x 10^° g/cc, Qc = 111.5, 

Qeq = 4.16, = 10^, Bmax = 7.22, giving rise to 

M = 1.95, re = 746, rp/re = 0.564 when KE/GE=0.1033, 
ME/GE=0.071; shown in Eig. [23b. 

(b) A white dwarf with pc — 2 x 10^^ g/cc, Qc = 40.56, 

= 4.25, = 10^, B max = 25.2, giving rise to 

M = 1.72, re = 428, rp/re = 0.689 when KE/GE=0.011, 

ME/GE=0.1133; shown in Eig. [23b. 

Eigures ESK and b, depicting the density isocontours 
of the above mentioned white dwarfs, indeed demonstrate 
the possible existence of white dwarfs with much smaller 
radii compared to their conventional counterparts. They 
could be, in fact, potential candidates for soft gamma-ray 
repeaters and anomalous X-ray pulsars with much smaller 
surface magnetic fields than that needed for a neutron star 
based model (magnetar). Now, one has to check the stabil¬ 
ity of such white dwarfs, in particular the one in Eig. 123b 
having higher pc, even though both models have acceptably 
small KE/GE and ME/GE. However, note that the thresh¬ 
old densities for the onset of several kinds of instability in 
general relativity in the presence of a high magnetic field 
would be expected to be different (probably higher) than 
those in Newtonian mechanics with a low magnetic field. 
Nevertheless, the value of F for the white dwarf in Fig. [23b 
need not be 4/3 throughout due to its higher Bmax- In spite 
of these cautionary notes, the above exploration argues for 
the existence of white dwarfs having re only an order or or¬ 
der and half magnitude higher than that of typical neutron 
star. Hence, they can be very good storage of a large amount 
of magnetic energy. 


8 SUMMARY AND CONCLUSIONS 

Magnetic field and rotation are important elements in the 
theory and phenomenology of white dwarfs, and realistic 
models must endeavour to include them. Differential rota¬ 
tion is expected to be a common feature in all kinds of star. 
As we have seen, it can significantly modify the structure of 
the star, and can cause the formation of concavities or hol¬ 
low regions near the poles, in special cases. The geometry 
of the magnetic field also plays a crucial role in determining 
the structure of the white dwarf, and toroidal and poloidal 
geometries introduce two entirely different kinds of defor¬ 
mation. We have computed sequences of equilibrium config¬ 
urations exhibiting magnetic fields of these two geometries, 
along with the inclusion of both uniform and differential ro¬ 
tations. We have adopted a GRMHD framework that has 
extensively been used by many authors to investigate neu¬ 
tron star structure, implemented numerically in the XNS 
code. We have modified XNS to work for white dwarfs, and 
the use of this code allows a better treatment of the non¬ 
linear Einstein equations in the presence of strong magnetic 
fields than in previous perturbative or Newtonian calcula¬ 
tions dm et al. ]H). 

For a white dwarf of fixed central density ~ 2 x 10^° 
g/cc, and magnetic fields of magnitude < 5 x lO^^G, we 
summarise below the key results from our parameter space 
study. 

• Differential rotation in the absence of magnetic fields 
can result in stars with masses of 1.4 — 1.8 M©, having equa¬ 
torial radii in the range 1200—1450km. The mass, equatorial 
radius, and the ratio of kinetic to gravitational energies in¬ 
crease with the increase of central (and, therefore, surface) 
angular velocity, apparently without an intrinsic bound. The 
largest surface equatorial angular velocities considered cor¬ 
respond to rotation periods of 1 — 2 seconds. There is a 
marginal decrease in the degree of central condensation, with 
mean density decreasing by about 30%. These stars show an 
oblate deformation, with the ratio of polar to equatorial ra¬ 
dius dipping as low as 0.6. In special cases, when the central 
angular velocity is greater than that of the corresponding 
uniformly rotating configuration at the mass-shedding limit, 
concavities or hollow regions develop around the poles. Slow 
variation of angular velocity in the radial direction, as com¬ 
pared to steep variation, causes larger deformations (always 
oblate), but less increase in mass. However, such flat rotation 
profiles cause smaller polar hollows, with the hollows disap¬ 
pearing altogether when the rotation profile approaches rigid 
body (uniform) rotation. 

• Toroidal magnetic fields in the absence of rotation give 
rise to stars in the mass range 1.4 — 2.3 M©, and equatorial 
radius range 1200 — 2330 km. Mass, radius and magnetic 
to gravitational energies ratio increase with the increase of 
maximum value of the magnetic field, which we have al¬ 
lowed to be as large as 3.4 x lO^^G. The surface of these 
stars deviates very little from a spherical shape, as indi¬ 
cated by low values, < 1.08, of the ratio of polar to equa¬ 
torial radii. However, the matter distribution in the interior 
shows extreme prolate deformation. The mean density falls 
by almost a factor of 10, indicating bloated outer regions of 
low density. When differential rotation is also present, the 
mass range is enhanced to 1.4 — 3.1 M©, with radii increas- 
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ing to nearly 3300 km. The corresponding magnetic field has 
maximum magnitude of 3.6 x 10^^G and the correspond¬ 
ing surface equatorial rotation period is about 10 seconds. 
Mass and radius both increase when the maximum held or 
the central angular velocity is increased, apparently without 
bound. However, peaks can be seen in the ratio of kinetic 
to gravitational energy and the surface equatorial angular 
velocity as functions of maximum held and central rotation 
rate respectively, indicating the possible onset of some kind 
of instability. The interior regions continue to show a prolate 
deformation, while the exterior tends to be oblate. When 
the differential rotation also causes polar hollows, the ratio 
of polar to equatorial radii falls to about 0.6. Differential 
rotation induced polar hollows are aggravated by toroidal 
magnetic helds. 

• Poloidal magnetic helds in the absence of rotation give 
rise to stars in the mass range 1.4 — 1.6 M©, and equatorial 
radius range 1200—1000 km. Mass and the ratio of magnetic 
to gravitational energies increase with the increase of maxi¬ 
mum value of the magnetic held, while the radius decreases. 
We have generally taken helds with maximum magnitude as 
high as 3.9 x 10^^G. Both the surface and the matter distri¬ 
bution show oblate deformation, with values of the ratio of 
polar to equatorial radii in the range 1 — 0.79. The mean den¬ 
sity increases by nearly 80%, indicating highly centrally con¬ 
densed and compact structures. When differential rotation 
is also present, the mass range is enhanced to 1.4 — 2.1 M©, 
and the equatorial radius range to 1400 — 1000 km. The 
corresponding magnetic held has maximum magnitude of 
5 X lO^^G, and the corresponding surface equatorial rota¬ 
tion period is about 2 seconds. There is an extreme oblate 
deformation, with the ratio of polar to equatorial radius de¬ 
creasing to about 0.55. The maximum value of the magnetic 
held generally increases with the increase in rotation rate, 
when the parameters determining the held are held hxed. 
Mass increases with both the maximum held and the an¬ 
gular velocity, apparently without bound. However, the ra¬ 
dius can either decrease or increase, depending on whether 
the held or the rotation effects dominate. Although both 
poloidal helds and rotation act in the same way with regard 
to deformation, namely by causing oblateness, the rotation 
causes increased size while poloidal held causes decreased 
size. If polar hollows are present, increasing the magnitude 
of the held results in a smoothening out of the polar hollows, 
due to the possible cancellation of centrifugal pressure and 
magnetic pressure. Furthermore, by increasing the central 
density and maximum magnetic held, super-Ghandrasekhar 
white dwarfs could be shown to be much smaller with the 
equatorial radius 400 — 750 km. 

The points noted above lead us to conclude that white 
dwarfs with highly super-Ghandrasekhar masses, i.e. M > 
Mch — 1.44 M©, are possible when rotation and magnetic 
helds are taken into consideration. 

Anticipating the existence of an upper mass 
limit for magnetised and rotating white dwarfs 
(|Das Mukhopadhyay I l2013al b we propose that the 
equilibrium solutions we have constructed here could be 
candidates for the progenitors of over-luminous, peculiar 
type la supernovae. The explanation of such supernovae 
requires highly super-Ghandrasekhar white dwarfs with 
masses in the range of 2.1 — 2.8 M©, and many of the 


conhgurations we have constructed in this work are within 
this mass range. 
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